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APPROXIMATE ANALYTICAL SOLUTIONS FOR HYPERSONIC FLOW OYER 

SLENDER POWER LAW BODIES 

By Harold Mirels 


SUMMARY 

Approximate analytical volutions are presented 
for two-dimensional (<r~0) and axisymmetric (a=l) 
hypersonic flow over blunt-nosed slender bodies whose 
shapes follow a power law variation. In particular , 
the body shape is given by r b — x m where r b is the 
transverse body ordinate , x is the sir earn wise distance 
from the nose t and m is a constant in the range 
2/ (<r + 5)<C in<C 1 . Both zero-order (M— > <» ) solu- 
tions and first-order ( small but nonvanishing values 
of lj( MS) 2 ) solutions are presented , where M is the 
freest ream Mach number and 6 is a characteristic 
body or streamline slope . The zero-order shock shape 
R 0 is similar to the body shape for these flows. The 
solutions are found within the framework of 
hypersonic-slcnder-body theory. 

The limiting case m —1 corresponds to a wedge 
(cr-- 0) or cone {<x=T) flow . The limiting case 
m— 2/(<r~\~3) corresponds to a constant-energy flow 
(r h =0, R 0 ~x 2/(<H_3) ). The latter cases are included 
so that the present study may be applied to all flows 
wherein the zero-order shock shape is given by R 0 ~ 
x m with m in the range #/(<r+S) <m<l . Flow 
fields associated with shock shapes having values of m 
outside this range are also discussed. For all values 
of in, except in -~1 } certain portions of the flo w field 
violate the hypersonic-slender-body approximations } 
i while other portions are consistent with these approx- 
imations . For m = 1, all portions of the flow field 
are consistent with the approximations. 

The approximate solutions are found as follows . 
The asymptotic form of the flow in the vicinity of the 
body surface is used as a guide to write approximate 
expressions for the dependent variables. These 
expressions exactly satisfy the continuity and 
energy equations and contain arbitrary constants 
which are evaluated so as to satisfy boundary condi- 
tions at the shock. The approximate solutions do 


not satisfy the lateral momentum equation except at 
the shock and ( for the first-order problem) at the body 
surface. 

The results of the approximate solutions are 
compared with numerical integrations of the equa- 
tions of motion for various values of m and y ( ratio 
of specific heats). Good agreement is noted , par- 
ticularly when m and y are both near one . The 
shock is relatively close to the body for the latter 
cases. Sufficient results are presented to evaluate 
the accuracy of the approximate method for various 
values of m and y. 

INTRODUCTION 

The steady-state equations of motion for hyper- 
sonic flow over slender bodies ran be reduced to 
simpler form by incorporating the “hypersonic- 
slcnder-body approximations” (c.g., refs. 1 and 2). 
The reduced equations are valid provided <$ 2 <C<C1 
and l/i\/5 <0(1), where M is the frcc-stream Mach 
number and <5 is a characteristic body or stream- 
line slope. Reference 1 has shown that, if the 
nondimcnsional streamwise coordinate is con- 
sidered as a nondimcnsional time, these reduced 
equations are identical with the full (exact) equa- 
tions for a corresponding unsteady flow in one less 
space variable. Forebody drag on a hypersonic 
slender body is equivalent to the net energy per- 
turbation (from the undisturbed state) ill the 
corresponding unsteady flow. 

References 3 and 4 have treated the constant- 
energy flow field behind the spherical “blast” 
wave which is generated when a finite amount of 
energy is released instantaneously at a point. The 
analysis assumes a very strong wave and is valid 
(for a perfect gas) until the decay of shock strength 
is sufficient to violate the strong shock assump- 
tions. The problem of planar and cylindrical 
blast waves was treated by a unified analysis in 

l 
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reference 5. In addition, the flow-field modifica- 
tions associated with more moderate shock 
strengths were found by a perturbation analysis 
(refs. 5 and fi) . The solution for the cylindrical 
blast wave was obtained, independently, in refer- 
ence 7. "Reference 8 has pointed out that, within 
the framework of hypersonie-slender-body theory, 
the hypersonic flow over a blunt-nosed flat plate 
(or circular cylinder) may be considered as the 
steady-state analog of the constant -energy planar 
or cylindrical blast-wave problems, respectively. 
The nose drag in the steady problem is equivalent 
to the finite energy which is instantaneously 
released in the blast-wave problem. The steady 
flow is not correct near the nose (where t lie 
hypersonie-slender-body approximation 5 2 <^<H is 
violated) and far downstream of the nose (where 
t lie approximation I/A/5<0(I) is violated as is 
the strong wave assumptions of blast-wave 
theory). However, useful results arc obtained 
for the intermediate regions (ref. 8). 

The blast-wave problems all exhibit flow 
similarity. That is, the flow fields at different 
times are similar, except for a scale factor on 
botli the dependent and independent variables. 
References 9 and 10 lmve observed that such 
similarity exists whenever the shock shape follows 
a power law variation (with sire am wise distance) 
provided the hypersonie-slender-body equations 
are considered in the limit as l/(A/5) 2 - >0. This 
led to numerical solutions of the hypersonic flow 
over slender blunt-nosed bodies. The effect of 
nonvanishing values of l/(il/5) 2 was also found in 
reference 10 by a numerical perturbation analysis. 

In the present report, approximate analytical 
solutions are obtained for both the zero-order 
(1 /(A/5) 2 — >0) and first-order (small but nonvanish- 
ing values of 1 /(A/5) 2 ) hypersonic flow over hi lint- 
nosed slender bodies. The zero-order solutions 
represent generalizations of the approximate ana- 
lytical solutions of the blast-wave problem which 
are presented in references 3, 4, and 5. The 
shock locations and pressure distributions indi- 
cated by the approximate solutions are compared 
with the values which result from a numerical 
integration of the equations of motion. Finally, 
some general properties of the hypersonic flow 
fields associated with power law shock waves are 
discussed. 

ANALYSIS 

The equations for hypersonic flow over slender 
bodies are summarized herein. These are then 



Figure I. — Physical quantities for study of hypersonic 
flow over blunt-nosed bodies. 7?o^Po(u) is shock shape 
in the limit as 1 /(JIT’S) 2 — 41. 


specialized to obtain the zero-order and first- 
order equations for hypersonic flow over those 
slender blunt-nosed bodies whose zero-order shock 
shape follows a power law variation. Finally, 
approximate analytical solutions of these equa- 
tions are obtained. 

HYPERSONIC-SLENDER-BODY THEORY 

The equations of motion for hypersonic flow 
over slender bodies (c.g., ref. 2) are summarized 
in the present section. 

Physical dependent and independent variables 
are barred herein (u f j, 7, etc.). Symbols are de- 
fined in appendix A. Figure 1 shows some of 
these quantities. Let 5 represent a characteristic 
body or streamline slope and L represent a char- 
act eristic streamwise length. Two-dimensional and 
axisymmetric flows are considered, with (x,r) 
and (5,7) being the streamwise and transverse 
coordinates and velocities, respectively. T o obtain 
the hypersonie-slender-body equations of motion, 
the following nondimensional quantities are intro- 
duced (following ref. 2): 

jr—7/L u= (5 — u J/u w 6 2 P^p/yMWp* 

r—rjLd v=v/u w S p=p/p« 

The body shape and shock shape are denoted by r b ~ 
F ft (r) and /?—/?(?), respectively, so that 

r 0 =h!L5 It=It/U (2) 

If these quantities are introduced into the equa- 
tions of motion and terms of order 5 2 are neglected 
(compared with one), the hypersonic-slender- 
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body equations are obtained. These are (ref. 2): 
Continuity: 


^ + ^ +<r pi’ =0 

dr dr r 


(3 a) 


r— Momentum: 


Energy: 



d(p/p y ) , d( 7 >/p0 „ 

dr “ r dr 


(3b) 

(3c) 


The boundary conditions are: 
At body surface : 



Upstream of shock: 

u m = v co = 0 
p m — 1 yd Pb 2 

P* — 1 


Downstream side of shod 


'7+1 


Vs 


_2 
7 + 


£{>-['/(£“)■]} 

! { L/ 1 . (Y- )]} 


_7+1 

7 “ 


(4a) 


(4b) 

(4c) 

(4d) 

(4e) 

(4f) 

(4g) 


Here <r=^0,l for two-dimensional and axisyminetric 
flows, respectively. This system of equations 
can be solved independent of the ^-momentum 
equation, and therefore the latter is neglected 
herein. The system of equations is consistent 
provided l/il/S<0(l). (Note that the right sides 
of eqs. (4c), (4e), and (4f) become infinite if Mb 
is permitted to go to zero as 5 2 goes to zero.) 
Tims, the conditions 1/d/S <0(1) and 6 2 <^<bL 
must be satisfied for the hypersonie-slender-body 
equations to be valid. 

FLOW FIELDS HAVING ZERO-ORDER SHOCK SHAPES 
FOLLOWING POWER LAW VARIATION 

The equations for hypersonic flow over slender 
bodies are specialized herein. The resulting 
equations give the zero-order and first-order 
hypersonic (low over those blunt-nosed bodies 
whoso zero-order shock shape follows a power law 
variation. These equations were previously de- 
rived in reference 10. 

Let /^(I) denote the shock shape for a given 
body in the limit as 1 / (A/ 5) 2 — >0 (i.o., zero-order 
solution). References 9 and 10 have shown that 
when /i*o(7) the flow fields are similar at each 


stream wise station (for 1/ ( A 2 — >0) . In general, 

the body shape is similar to the shock shape for 
such flows, in which case r+7)^7 m . Values of m 
in the range 2/(<r+3)</m<l correspond to bodies 
having an infinite positive slope at the nose. 
The limiting value m— 1 corresponds to flow over 
a wedge (cr = 0) or cone (<r~l), while m = 2/(cr + 3) 
results in a body shape r 6 (/)— 0. The latter flow 
may be interpreted as that over a blunt-nosed 
fiat plate (<r = 0) or blunt-nosed circular cylinder 
(<r— 1), (As in refs. 9 and 10, the m — 2/(<r + 3) 
case is referred to as the “const ant -energy case.”) 
The present report is primarily concerned with 
values of m in the range 2/(<r+3) < l. 1 These 

flows violate the hypersonie-slender-body assump- 
tion 5 2 <( <C 1 11 1 7 — 0 (except for m = 1) so that the 
resulting solutions are not expected to be valid in 
this region. _ 

For the remainder of this report L is taken to 
be the slreamwise length of a given body and 5 is 
defined to be 

5=11,(1)11 (5) 

Then, the nondimcnsional zero-order shock shape 
is given by 

I?o=r m (0) 


An alternate characteristic slope based on body 
thickness at x L is 


^6 = 


+ +) 

L 


^5r b (\)=b 


r b (x) 
H oW 


( 7 ) 


These quantities are indicated in figure 1. For 
problems where the body shape is initially speci- 
fied, 5 b is known immediately while b is found as a 
consequence of the solution. 

New independent variable's are now introduced 
according to the relations 


1 

T l =rJH l) =r/j m J 

so that 

c) c) mij 2) ] 

< 57 — T&n 



(8) 


( 9 ) 


i Values of m outside this range give rise to flow fields, portions of whicli are 
physically realistic and consistent with the assumptions of hypersonie- 
slender-body theory. These flows are discussed in the section entitled 
GENERAL CHARACTERISTICS OF FLOW FIELDS ASSOCIATED 
WITH POWER LAW SHOCKS. 



4 


TECHNICAL REPORT R-lo— NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


At the zero-order shock location 77 — 1 , and at the 
body 77 is also constant and is denoted by 7 ^. 
Note also, 5 & = 5? Following references 5, 6 , 
and 10 , a small perturbation parameter € is intro- 
duced : 

e=(?-™JM8y 

d e/dr —2(1 — in') ef£ 

The boundary conditions at the shock suggest the 
following forms for the dependent variables: 


v — to £ m 1 (y>o + t<p\ ) 

(11a) 

p = m 2^m-l) (Fo+tFi) 

(lib) 

p=4 / o+(4 / i 

(He) 

7?=P (l+ea.) 

(lid) 


where <p, F, and \j/ are functions of tj only and ciy 
is a constant wliicli is initially unknown. 

An alternate shock shape parameter 0 can be 
defined according to the relation 



The range 2/(o-+3) < m < 1 then corresponds to 
the range 1 >$> 0 . Note, 0—1 corresponds to the 
constant-energy case. 

Substituting these quantities into the equations 
of motion and collecting terms of order c° and e 1 
then define the zero-order and first-order approxi- 
mations, respectively, for the hypersonic flow over 
a power law body. The choice of the variables is 
such that the zero-order and first-order equations 
are functions of 77 only. These results, which were 
previously obtained in reference 10 , are summa" 
rized as follows. 


Zero-order approximation. Equations (3) yield: 
Con ti unity: 

(<Pu — ^'oV ? o“h <r — 0 

V 

(13a) 

17 — Momentum : 


(<p Q — 7 ?) + ^ — p) 0 

(13b) 

Energy: 

( 1 3c) 


where primes indicate differentiation with respect 



to 77 . The boundary conditions at tj=1 are 


*,(1)=F 0 (1)=2/(<H-1) 
^o(1) = (y+1)/(7— 1) 


( 14 ) 


Equations (13) and (14) define the flow field com- 
pletcly. The body location is defined by the 
tangency condition (eq. (4a)), which becomes 
(since r b = 7j b Tf 0 ) 

^o(vb) = Vb (15) 

Derivatives of the dependent variables, at r\= 1, 
are given in appendix B (for later use). 

First-order approximation. — Equations (3) yield, 
respectively, 




V>i 


+ 


k («^+l)fS — ( spo — v) 7^ 7^=0 (10a) 

L r o J ro 


-3 ^=o 


<P \ 


(1Gb) 


Vi+fe*. - n ) ^+(^+7) Fp° 

(16c) 


The boundary conditions at 77 -= 1 (see appendix 
C) are 



_r -2 ( l+ c+\ y-| 


i-«i 0) 

t_7 — 1 V 2 Vj 


— 

^(i)“— 1 [i + ((r+i)0]^ 

dl 

== 

615 bi^L\ 

(17a) 

*.o) r 
^«(i) L 







= bi&— & 24«1 

(17b) 

L(i) r 
F 0 ( 1) L 

^r’O+T*) 8 ] 



— ' ^o(l) - 2[1 + (c+ 1)0] ^ a x 


=b^~b u a L (17c) 
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The tangency condition at the body has already 
been satisfied (eq. (15)) so that the boundary 
condition on <p x at the body surface is 

<Pi = 0 O 8 ) 

For a given body, a x must be determined such 
that equations (16) and (17) yield a solution 
satisfying equation (18). To avoid trial-ai id-error 
choices for a { in a numerical integration of equa- 
tions (16) and (17), it is advisable to decompose 
the dependent variables and boundary conditions 
into two parts, one independent of a x and the 
other proportional to a x . That is, each dependent 
variable is expressed as 

( ) x =( ) M + ( ) lt2 a, (19) 

For example, mid so forth. The 

solution for ( )i,i and ( )i, 2 can be obtained 

independent of a X) and the final solution is given 
l) 3 r equation (19) with a x found from 

a x = — 

Such a procedure was used in references C and 10 
and is permitted because of the linearity of 
equations (16) and (17). This procedure was 
also used herein to get- additional numerical 
integrations of equations (16) and (17). When 
getting approximate analytical solutions of these 
equations, it is possible to satisfy equation (18) 
without resorting to equation (19). 

EXPRESSIONS FOR SURFACE PRESSURE, SHOCK SHAPE, AND 
DRAG 

Before continuing with the solution of the zero- 
order and first-order problems, it is useful to 
develop expressions for surface pressure, shock 
shape, and drag for these flows. 

Within the framework of hypersonic-slender- 
body theory, the local pressure coefficient is given 
by (from eqs. (1)) 

r = —op P j) — — - — "] 

” _ i- - 2 l p 7(WJ (20) 

2 

The local surface pressure coefficient for the zero- 
aud first-order problems is then given by either of 
the following expressions (from eqs, (lib) and 
(20)): 

= 2KM f , , |T iM 

(drjdx ) 2 nl \ 

L 1 "I ( 

7m 2 / ( ’ 0 ('i?,l J (A/<5 4 ) 2 \ 



or 


(dn u /d7y 


■■2F 0 ( Vl ) 


1 

ym 2 F 0 (ri 


( l+ K 

5 ] 




■ 1 / 

7yc-™) 

. WV 

u J 


(21b) 


The equivalence of equations (21a) and (21b) can 
be seen by noting T b =r] b B [X . Equation (21b) is 
particularly useful for the jS= l (constant energy) 
case since r b — 77*= 0 therein. 

The zero-order and first-order shock shapes are 
given by (from eq. (lid)) 


J?Cr) , 

r b (x) Vh 




l+fl-i 


n» 

(did,) 


or 


JM 

/?„(/) 


=i+ar 




jr\ 2 (l-m) 


(0 


(22a) 


(22b) 


Again, equation (22b) is particularly useful for 
the 0—1 (*ase, as is discussed later in this section. 

Numerical values of the quantities 1 ?&, (^?^) > 

Fi(ri b ) f and a u which appear in equations (21) and 
(22), are tabulated in tables I to I\ for various 
values of < 7 , 7 , and 13. 

The forebody drag can be found by integrating 
the surface pressures. Tf D(jr) is the forebody drag 
up to station / and g = Pw^i/2, t lien the appropri- 
ate inf.eeral is 


^ = 2.° P** 

<Z Jo 


C,. »(?*)' dr„ 


(23) 


Noting r 6 =S i ,U'", substituting equation (21 u) into 
equation (23), and integrating yield for 
m>2/(<r+3) 

D(x) __ = 2F u (v„ ) f r (g+3) ~ 2 

2ir a m 3 q8l{LS 6 Y + ' nl \ m(<r+3) — 2 

, r f * +l) m*) i , _l , f \ ( 24 ) 

1) 7'» 2 ^o(>?o) J (Alia)* f 

The over-all forebody drag coefficient, referenced 
to the cross-sectional area of the base, is then 


„ _ D(L) 2 / 1 

Cw= 2^[n(r)r +i “[«(^+3)-2i ni \ 

, m(<r f 3)-2 I - Fi(r}fi) 1 I nl \ (2S) 

w(<r+l) LF 0 (i».) 7W z F 0 (j7 ft )J (Mdo) 2 j 
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-APPROXIMATE SOLUTION OF ZERO-ORDER PROBLEM FOR AND COMPARISON WITH 
RESULTS OF NUMERICAL INTEGRATION OF EQUATIONS OF MOTION 


y 

0 

Approximate solution 

Numerical 

integration 

r„ 

A. 

Vb 

Foivb) 

Vb 

FuO n) 

1. 15 

0 

0 

1. 255 

0. 930 

0. 930 

0. 930 

0. 930 


/3 

2. 93 

I. 535 

. 891 

. 766 

. 893 

. 761 


/2 

5. 18 

L 800 

, 852 

. 672 

. 852 

. 675 


h.L 

7. 16 

2. 13 

, 801 

. 598 

. 803 

. 611 


Vi 

8. 52 

2. 68 

. 710 

. 520 

. 716 

. 546 



6. 36 

3. 80 

. 513 

. 442 

. 535 

. 483 


1 

. 465 

7. 67 

0 

. 412 

0 

. 415 

1. 4 

0 

0 

L 556 

0. 833 

0. 833 

0. 833 

0. 833 


/3 j 

H 

. 962 

1. 078 

. 760 

. 679 

. 759 

. 666 


1. 240 

1. 779 

. 695 

. 584 

. 695 

. 581 


% 

L 296 

1. 888 

. 619 

. 504 

. 623 

. 518 


% 

L 167 

2. 04 

. 499 

. 415 

. 513 

. 454 


% 

. 817 

2. 30 

. 284 

. 320 

. 333 

. 390 


i 

. 417 

3. 50 

0 

. 316 

0 

. 325 

L 67 

0 

0 

1. 788 

0. 749 

0. 749 

0. 749 

0. 749 


X 

. 535 

1. 788 

. 660 

. 605 

. 658 

. 587 


H 

. 644 

1. 770 

. 586 

. 532 

. 585 

. 507 


X 

. 655 

1. 762 

. 505 

. 432 

. 509 

. 146 


% 

. 608 

1. 752 

. 385 

. 340 

. 404 

. 386 


% 

. 504 

3. 738 

. 186 

. 239 

. 248 

. 326 


i 

. 375 

2. 49 

0 

. 250 

0 

. 264 


APPROXIMATE SOLUTION OF ZERO-ORDER PROBLEM FOR a =\ AND COMPARISON WITH 
RESULTS OF NUMERICAL INTEGRATIONS OF REFERENCE 10 


y 

0 

Approximate solution 

Numerical integra- 
tion (ref. 10) 

c„ 

Do 

i7b 

FoOh) 

Vb 

Fo(vb) 

1. 15 

0 

- 2. 53 

2. 05 

0. 965 

0. 947 

0. 965 

0. 94S 


X 

-1-1. 636 

1. .317 

. 945 

. 774 

. 945 

. 775 



3. 64 

1. 042 

. 924 

. 682 

. 924 

. 688 


% 

5. 53 

1. 981 

. 897 

. 60S 

» 898 

621 


% 

7. 1 1 

2. 53 

. 840 

. 533 

. 845 

. 553 


% 

5. 68 

3. 62 

. 724 

. 453 

* 735 

* 484 


i 

. 465 

7. 13 

0 

. 410 

0 

.411 

L 4 

0 

-1.019 

2. 11 

0. 915 

0. 872 

0. 915 

0. 875 


}i 

+ . 383 

1. 228 

.875 . 

. 704 

. 875 

. 696 


)i 

. 757 

L 485 

. 839 

. 611 

. 839 

. 607 


% 

. 930 

3. 631 

. 795 

. 529 

* 796 

» 538 


% 

. 949 1 

1. 798 

. 719 

. 438 

. 725 

. 467 


% 

. 757 

2. 03 

. 561 

. 337 

» 589 

«. 392 


1 

. 417 

2. 92 

0 

. 302 

0 

. 311 

1. 67 

0 

-0. 625 

2. 16 

0. 870 

0. 805 

0. 870 

0. 81 1 


X 

+ . 170 

3. 161 

. SI 0 

. G45 

. 819 

. 634 


H 

. 365 

1. 405 

. 776 

. 551 

. 776 

. 544 


% 

. 452 

1. 450 

. 726 

. 469 

* 727 

* 474 


% 

. 483 

1. 457 

. 64 4 

. 370 

. 652 

. 403 


% 

. 459 

3. 436 

. 480 

. 257 

*. 518 

* 326 


l 

. 375 

1. 867 

0 

. 255 

0 

. 241 


- Numerical integrations for 0=H, J6 were not given in ref. 10 and were found as part of the present study. 
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TABLE III— APPROXIMATE SOLUTION OF FIRST-ORDER PROBLEM FOR <r=0 AND COMPARISON 

WITH RESULTS OF NUMERICAL INTEGRATION 


y 

0 

Approximate solution 

Numerical 

integration 

At 

/?> 

c, 

(h 

F7.ru,) 

d\ 

F ’.(») 

1. 15 

0 

1. 935 

0 

-13. 33 

1. 000 

1. 800 

1. 000 

I. 800 


X 

4. 47 

0516 

-27. 6 

1. 342 

2. 78 

1. 34 

2. 79 


X 

6. 91 

+. 0550 

-35. 4 

1. 453 

3. 03 

1. 43 

3. 16 



9 20 

. 226 

-41. 8 

1. 474 

2. 95 




/s 

% 

1 1. 76 

. 385 

-48. 3 

1. 415 

2. 57 

1. 31 

2. 77 


% 

15. 11 

0457 

— 54. 3 

1. 194 

1. 999 




1 

14. 82 

+ 6. 85 

-60. 3 

. 963 

. 923 

1. 03 

. 910 

1. 4 

0 

1. 857 

0 

-5. 00 

1. 000 

I. 548 

1. 000 

I. 548 


X 

3. 06 

1035 

-11. 28 

1. 213 

2. 18 

I. 21 

2. 17 


X 

4. 24 

-. 1135 

-14. 15 

1. 230 

2. 22 

1. 21 

2. 25 



5 14 

1110 

- 16. 18 

1. 192 

2. 06 




7° 

% 

5. 98 

-. 1430 

-18. 21 

1. 113 

1, 764 

1. 07 

1. 78 


7 A 

6. 82 

— . 386 

-20. 0 

. 975 

1. 397 




/a 

1 

6. 01 

2. 94 

-22. 7 

. 992 

. 781 

. 965 

. 799 

1. 67 

0 

1. 799 

0 

-2. 99 

1. 000 

1. 348 

I. 000 

1. 348 


X 

2. 27 

-. 1000 

-6. 99 

1. 140 

1. 778 

I. 14 

1. 77 


X 

2. 95 

— . 1237 

-8. 55 

1. 130 

1. 750 

1. 11 

1. 78 


Vs 

3. 43 

1323 

-9. 70 

1. 092 

1. 599 




X 

3. 85 

1644 

-10. 83 

1. 035 

1. 367 

. 991 

1. 39 


% 

4. 22 

-.309 ! 

-11.90 

. 956 

1. 092 




/o 

1 

3. 96 

+ 1. 813 

-13. 22 

. 940 

. 700 

. 930 

. 707 


TABLE IV— APPROXIMATE SOLUTION OF FIRST-ORDER PROBLEM FOR <r=l AND COMPARISON 

WITH RESULTS OF NUMERICAL INTEGRATION 


y 

0 

Approximate solution 

Numerical integration 


Present results 

Kef. 10 

At 

B x 

c, 

a A 

F,(vb) 

flj 

Fi (vi,) 

fli 

F,M 

1. 15 

0 

I. 163 

0. 00819 

-13. 16 

0. 489 

1. 10 

0. 489 

1. 10 

0. 455 

1. 10 


x 

4. 03 

. 229 

-34. 1 

. 883 

2. 55 

. 885 

2. 50 

. 829 

2. 43 


u 

7. 33 

. 696 

-48. 7 

1. 083 

3. 28 

1. 08 

3. 37 

. 982 

3. 05 


X 

10. 78 

1. 545 

-61. 7 

1. 204 

3. 52 

1. 17 

3. 84 




X 

14. 53 

3. 25 

- 76. 4 

1. 288 

3. 21 

1. 18 

3. 79 

1. 15 

3. 53 


% 

17. 00 

6. 61 

-93. 2 

1. 377 

2. 26 

1. 10 

2. 93 




1 

17. 25 

3. 50 

-77. 5 

1. 449 

1. 064 

1. 23 

1. 32 

1. 07 

1. 35 

1. 4 

0 

1. 027 

0. 00909 

-4. 83 

0. 477 

0. 899 

0. 476 

0. 918 

0. 396 

0. 92 


X 

2. 74 

. 1231 

- 13. 39 

. 802 

2. 03 

. 807 

1. 97 

. 077 

1. 93 


X 

4. 57 

. 315 

— 18. 83 

. 931 

2. 50 

. 932 

2. 49 

. 793 

2. 38 


% 

6. 21 

. 676 

-23. 4 

. 991 

2. 59 

. 976 

2. 67 




% 

7. 77 

1. 391 

-28. 5 

1. 027 

2. 36 

. 976 

2. 51 

. 870 

2. 36 


% 

8. 73 

2. 81 

-34. 2 

1. 075 

1. 796 

. 964 

1. 95 




1 

7. 67 

1. 372 

-28. 3 

1. 231 

. 954 

. 992 

1. 14 

. 937 

l"07 

1. 67 

0 

0. 968 

0. 01858 

-2. 84 

0. 467 

! 0. 785 

0. 465 

0. 783 

0. 350 

0. 78 


/3 

1. 996 

. 0966 

-8. 11 

. 754 

j 1. 692 

. 762 

1. 63 

. 575 

1. 61 


X 

3. 19 

. 225 

-11. 26 

. 856 

2. 03 

. 863 

2. 00 

. 603 

1. 93 


5 A 

4. 21 

. 472 

- 13. 87 

. 903 

2. 09 

. 900 

2. 11 




% 

5. 18 

. 952 

-16. 69 

' . 933 

1.919 

. 911 

1. 98 

. 732 

1.76 


% 

5. 83 

1. 906 

-19. 85 

. 974 

1. 606 

. 922 

1. 57 




1 

5. 10 

1. 056 

-17. 16 

1. 118 

. 823 

. 969 

. 931 

. 820 

790 


For 7 = 1.405, 0=0, <r = l, tltc results of ref. 2 indicate a, =0.47 and Fi (t?&) =0.91 (found from cone results for M6b=<*>, 3.988 therein). 
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For the constant-energy case, 0 — 1 or m— 
2 /(<r+ 3 ), the integration of equation (23) is invalid. 
However, the drag for the latter case can be ex- 
pressed in terms of a momentum contour integra- 
tion. A momentum contour integration indicates 
that, the forebody drag up to any station x equals 
the net energy perturbation of the transverse flow 
(per unit x) at that station (e.g., ref. 10 ). The 
energy perturbation is taken to be the departure 
from the free-stream value. Thus, if E is the 
energy perturbation per unit mass at any point, 


s.,(T-T J+ !- 7 4t(B;)+? 


then 


B(J)=2r 


J 'rG) 

rtCr) 


pEr° dr 


where 


P E ^~ I (^.+- 2 1 


+ 


1 

(A/5 ) 2 _ 


7—I 


ym* 


0 


(& + Sr)]} 


(26) 

(27) 


(28) 


Noting r = 5 /.£ m 77 , d r — 5L£ m dq and integrating be- 
tween the limits r\ b to 1 + etfi yield 


D(x) 


4 7r tr m 2 g(5/ v ) < 7 + 1 5 2 

gm (<r + 1 ) 


■m (ff+ 3 ) — 2 


f(Z ?+ i , dj? 


f4a i 


t (t-1) (Msy\y+\ ym*{<T+\) 

+ ^ dT? } (29) 


wherein use was made of the relation 



l/(cr-fl). Equation (29) is applicable for all m 
and is thus more general than equation (24). 

Consider the zero-order drag approximation 
(first term on right side of eq. (29)). It is seen 
that the dependence on £ disappears for m— 
2 /(<r+ 3 ) so that the drag D(x) and, therefore, the 
transverse energy jump discon tinuously at x ~0 
and are constant for x)> 0 ; hence, the origin of the 
term constant-energy case for this value of m. 
If the first-order approximation is to be included, 
equation (29) shows that in order to have a con- 


2 By employing the zero-order stream function (see next section), 
0 q/(<H- 1 )»j*, so thatj^ twf d 77 =ri— 0 o(tjb)]/(H-l). For body shapes de- 
fined by 0 o(j?fc) = O, the integral becomes l/(<r+l). 


s tan t-en orgy flow the value of c?i must be such as to 
make the coefficient of £ mC<r+1) identically zero 
in equation (29). Such constant-energy flows 
can be used to estimate the shock shape and 
pressure distribution on blunt-nosed flat plates or 
circular cylinders (see following sketches). The 
correspondence between these two flows (sketches 
(a) and (b)) breaks down at x — 0 because of the 
finite thickness of the plate or cylinder. How- 
ever, hypersonic-slender-body theory (which is 
the starting point of the present analysis) is 
inapplicable near j =0 anyhow. These two flows 
are expected to be in essential agreement in the 
intermediate regime wherein the values of x are 
neither too small nor too large to invalidate 
hypersonic-slender-body theory. 



(a) Constant -energy flow. Drag impulse at 7 = 0. 


r 



(b) Flow over blunt-nosed plate or circular cylinder. Drag 
impulse at x= 0. 

If the nose drag (i.e., drag impulse at x= 0 ) is 
known in a constant-energy problem, it is possible 
to express the shock sliapcand pressure distribution 
as a function of this known nose drag. The 



APPROXIMATE ANALYTICAL SOLUTIONS FOR 
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procedure is as follows. For a constant-energy 
flow, equation (29) becomes 

Dn 

4*°m*q[R 0 (j)]° + '[R 0 (*)F]* 

=£ (^5r+| Vito) -rdus I (30) 

whtTc T) n is the known nose drag and the substitu- 
tions L=x and d=F 0 (x)pc have been made (since 
the equation is independent of the choice for X). 
Let C Dw be the nose drag coefficient: 


reference 5 (table V(b)). The results are expressed 
as (( 7 + 1 ) (t 2 — 1 )/ since this is a slowly varying 
function of 7 and <r. Values of a x are listed in 
tables III and IV for < 7 — 0,1 and 7=1.15, 1.4, and 
1.67. 

The corresponding expression for surface pres- 
sure coefficient for the constant-energy case is 
(from eq. ( 21 b)) 


Oj?,b 

(d/? 0 /dx) 


^=2F o (0)Il+R^ 

,/dx ) 2 l L^o( 


(cr+3) 2 ' 
47^o (0). 


C n .=; 


By 

Vq(r N ) c + l 


where r N is the half thickness or radius of the flat 
plate or circular cylinder, respectively. The shock 
shape can then be expressed as (from eq. ( 22 b)) 

MJ^/i +a Jf r J-Js-T> (32a) 

TV TV \ L M It 0 (x) J / 


where (from the ratio of eqs. (30) and (31)) 

5o(j) 1 r (^+3)^ 

?V 2 L I J VV/ 


(32b ( 


Equations (32) give the shock shape Ii(x)jr N ts a 
function of x[r N with C Dn > /, and a x as parameters. 
For a given nose configuration, Cp N can be esti- 
mated by using methods described in reference 11 
and in the references noted therein. Numerical 
values of I (defined by eq. (30)), obtained from 
an integration of the zero-order equations, are 
listed in table V(a) for various 7 and cr. These 
are in good agreement with values reported in 


where I{ 0 (x)[r N is given by equation (32b). 

FORMULATION OF ZERO-ORDER PROBLEM IN TERMS OF 
STREAM FUNCTION 

The zero-order and first-order problems could 
have been formulated in terms of a single depend- 
ent variable, the stream function, instead of the 
three dependent variables p, p, and v. Such a 
formulation is convenient for obtaining asymp- 
totic solutions for the flow in the vicinity of the 
body. These asymptotic solutions are useful when 
numerically integrating the equations of motion 
and when approximate analytical solutions are de- 
veloped. Hence, the zero-order problem is formu- 
lated in terms of a stream function herein, and 
asymptotic solutions are obtained for 77 near 
The continuity equation (3a) is satisfied by a 
stream function \j/ defined such that 

jV B= _5|/ax ( 34 ) 

r" c )r tyjdr 

The energy equation (eq. (3c)) shows that j>lp y is 

constant along a streamline (except for the dis- 


TABLE V. — EVALUATION OF 7= f'( 


F» ,1 „ > 

7 +o ^orQ 

7 1 Z j 


• d, AND COMPARISON WITH REFERENCE 5 


(77= J 


[Tabulation is In terms of (a-[-l)(7 2 — DI-] 

(a) Present results, based on numerical integration of (b) \ alues based on Table III of reference 5 
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continuity at the shock). By considering flow downstream of the shock, a function co (of fy) can 
be defined according to the relation w(^) =p/p T so that 

P = <*P y (35) 

For prescribed shock shapes, the functional dependence of w on ^ can be determined explicitly by 
considering conditions at the shock. Substitution of the previous equations into the momentum equa- 
tion reduces the problem to that of determining the single dependent variable 

In the present section, only the zero-order problem is considered. An appropriate form for is 


[!/(*+ l)]r ( " +1) 0o 

where 6 0 = d 0 (r}) and 0 o (l) = l. At n=l, the boundary conditions on p 0 and p 0 give 

2 


w 0 = 


But, also at 1, 


- ( - — -Yw 2 £ 2(m “ 1) 

7+1 \7+l/ * 


^0 


a+l 


(36) 


(37a) 


(37b) 


Fliminating £ between equations (37a) and (37b) and substituting for according to equation (36) 
show that 

2 m 2 / y — 1 




7+1 \7+ 


\y \r fv+l X]- fi 


(38) 


Therefore, p,„ v 0 , and p„ can be expressed as functions of f and 0 O (from cqs. (34), (35), (36), and (38)). 
Substitution into the momentum equation (eq. (3b)) yields 


0 . <r 00 , (y+l) 1 ' +1 (ff+l)' y V 1+<r ' 1 ' 1 >0'+0 

r^-i 

2 ' 

(„ , *+1 


00 00 V V 00 27(7- 1 ) T (0o) y 

(' + 2 

0 00 ) 2 , ( 7 +P T+ 1 (<H-i) T V 1+ ^~% 

i +<, i 

(tr + 1)0 O 
_ V0o _ 


2y(y-l) y (ff') y | 


A V 
/ voL 


(39) 


with the boundary conditions 
<?o(l) = l 

>+l)( 7 +l) 


«£(!)= 


7—1 


As 0 O >0 o (tj s )=O (i.e., 77 — >tj 6 ) , the last term in 
the numerator and denominator of the right side 
of equation (39) can be neglected (at least for 
/3< l). 3 Equation (39) can then be written 


On the body, 0oO?»)=O. 

The dependent variables of the previous section 
are related to 0j by the relations 


#o_<r „ 

Oo 7 So V 


(41) 


<Po=V—(.<r+ 1) X7 
“0 

, 1 00 


cr+1 if 


(40a) 

(40b) 


Integrating gives 

Oo^Koiv'+'-vl+'y'^-V (42) 

where K 0 is a constant. Substitution of equation 
(42) into equation (39) verifies the neglect of the 
last terms in the numerator and denominator for 
77 near An improved asymptotic solution can 
be obtained by substituting equation (42) into 



F ^ W 

' 1 7-1 0oY 

(40c) 

these terms. 

This gives for /3<y 1 and 

P=h 


U ~7+l 

,v + 1 7+1 r) 

respectively (recalling ti b =0 for /?=1), 


Also, 




'"1 

b 

<0. 
co 1 

(<r-\~l)pK 0 7fl 



■Cj — 

1 1 
$ 5, 

00 00 
00 O’o 

(40d) 

( 7^o 1 i 

_2(y-P)F 0 ( Vlj )J 

Xv"(v' ,+l —nl +1 ) 0,{y ~ 0) 

0<1 


S-IS-, 

II 

,0o__v 
00 V 

(40e) 

— A 

7+tr 

f Ko[y{a +3) — 2] "1 7-1 

L27 2 (7-1)E 0 (0) P 

0=1 


3 

F 0 7 00 00 V 


( 401 ) 


3 The limit is associated with the a= 1 case (see discussion in section 
entitled GENERAL CHARACTERISTICS OF FLOW FIELDS AS- 
SOCIATED WITH POWER LAW SHOCKS). 


APPROXIMATE ANALYTICAL SOLUTIONS FOR HYPERSONIC FLOW 


OVER POWER LAW BODIES 1 1 


Integrating the expression for 0<1 and then 
substituting into equations (40a), (40b), and 

(40e) give, for 0<1, 


e^K 0 W+ l -r,l + r ,iy -^\_ 


n- 


fi ijJ-' 




2(2 7 -jS)F 0 W 

_ 7 ? » + i )7/(y-^ > ] ( 43a) 


V^O" 


'I'o- 


t-zs^+'-^+T, 

— — L 1 - 




ft 


•] 


2(2 7 -0)F o (i?») 0 

(43b) 


^ [v’ +l -nl + T (y -^ |~1- 

y—p L 


fill- 


(2y-p)F 0 (; 




(43c) 
(43d) 


Integration then gives 
0o~ Ao(» 7 " +i — vl +} ) [l 

, , , * w x - ^«(^ +, -^ +, ) 4 l (431) 

6 7 (v+l)AT 0 ^ Fo(»;(,) J 

■ [' e °J 

(43j) 
l] (43k) 

>?] (430 


<f>0 — V ' 


U' +1 -^ +] 


'[■ 


2(<r+l) 7 M'Fo(r, 4 ) 




ft o«F 0 (ij») [l- 


Terms of order 0“ are neglected in the brackets 
of the previous equat ions. Integrating the expres- 
sion for 0=1 and then substituting into equations 
(40a), (40b), and (40c) give, for 0=1, 


6u^K 0 t] 


<Po—V' 


0 + 1)7 

7 — 1 


((7 + l)X- 0 n 


2 7 +g~ 

7-1 


(43c) 


"2(41)^WW 

Terms of order 0o are neglected. 

These asymptotic solutions (eqs. (43)) are useful 
in numerically integrating equation (30) (or eqs. 
(13)). The boundary conditions at tj= 1 permit 
one to integrate from 77= 1 toward = How- 
ever, equation (39) is singular at r} b , and the 
numerical process breaks down. But, the asymp- 
totic solutions can he used to carry the solution 
to For example, if <p 0 is eliminated between 
equations (40a) and (43b), the resulting equation 


r 2 7 (2 7 +<r-l)F 0 (0)J 


order, is 



r i 

(43f) 

r, 

7 (<r+l) 0 a V j 

fir, l ~°0 o T 

L 1 ~VF 0 (0j v e °_ 

'7s~’7-S 1- 

y—fi y6'o L 

1 2(2y—fi)F 0 J 


1/ (<r+l) 


a + 1 


_yKo T-ir . 3 y + «^(7+1)- 1 
fo ~ 7 -l l v L + V(2r+v-l)Fo(0) ^ °J 

(43g) 


Similarly, equations (43d) and (43h) give 


F 0 ~F 


o(0) [l 


37+<ry— 2 


2y(2y+a-l)F c 


(43h) 


F 0 {vb) 

F /. fi v l ~%\ 
~H 1_ 2 F 0 ) 

fi<l 

W 1 1 

37+0-7—2 + 

0=1 

Tq I I 

2 7 (2 7 +<r-l) Fq J 


Terms of order (ij ,_<r 0o) 2 are neglected. Note, from 
equations (40c) and (42), that K 0 is related to 
F 0 (v>) by 


/(T-/S) 


*-(?+ «*>] 

For 0 = 0 the correction terms vanish in equa- 
tions (43a) to (43d). A higher order approxima- 
tion can be obtained, for 0=0, by writing equation 
(39) as 


i V '( V *+*— v 

e ' 0 v Lyvl*Fo( Vb )r {ri v 


ff +i) 


The asymptotic solutions are also useful as guides 
for setting up approximate analytical solutions of 
the zero-order equations. This is done in the 
next section. 

APPROXIMATE ANALYTICAL SOLUTION FOR ZERO-ORDER 
PROBLEM 

An approximate analytical solution for the zero- 
order problem can be obtained in the following 
manner. Equations (43b), (43f), and (43j) sug- 
gest that <p 0 may be approximated by the expres- 
sion 

-,=-(*=£) 1 or +, -^ +l )[i-C,Or H 

V 7 ~vl +1 ) D °] (44a) 


<Po 
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where C 0 and T) y) (as well as are as yet unknown 
constants. Differentiation of equation (44a) 
shows 

^+T 9=(<r+1) [&=* CoffVH) (ir M -*r H ) n d 

v L7 y J 

, (44b) 

<Po + a — /—y ( u( D o 

+ l)D 0 i;'(^ +1 -i ? f +1 )®o-i (44c) 

Substituting equations (44a) and (44b) into the 
continuity equation (eq. (13a)) permits the latter 
equation to be integrated. The result is 


Any form for <p () — r) which permits equation (47) 
to be integrated in closed form will permit equa- 
tions (13a) and (13c) to be integrated in closed 
form. Similarly, equations (45) and (46) could 
have been deduced from equations (40b), (40c), 
and (48). 

The constants C 0 , D 0j and vo will now be evalu- 
ated. The cases 0< 1 and 0=1 are treated 
separately. 

Case 0<C 1. — Define the following known quan- 
tities: 


p 1 y 

° 7+1 7-0 


(49a) 


i» [i-c,(i -gC 1 )” 0 ] 

*° (1) (l 

X- y (45) 

Similarly, integration of the energy equation (eq. 
(13c)) yields 

_.7(7>o+1) 

i-Cod-qC 1 ) 00 1 ,, n 

The constants C 0j D 0 , and 77 6 will be determined 
so as to satisfy the boundary conditions on ^ 0 ( 1), 
V?o(l), and <^'(1) given by equations (14), (Dl), 
and (B4). The resulting solution for <p Q , F 0 , and 
will have the correct values of these func- 
tions and their first and second derivatives at 
V— 1 and will satisfy the auxiliary' condition 
<Po{Vb) In addition, the continuity and energy 

equations are identically satisfied, hut the momen- 
tum equation is not (except at t)— 1). 

Note that t lie expression for associated with 
equation (44a) is (from eq. (40a)) 


0 0 =erp 


+ (*+i) f ----- 

L Jv <Po—yJ 


. [ l-r„d- 

(i-„e‘+ 

K (Y +1 -^ +1 +* 

[l — C«(tj' +1 — 


So— ( 0 ._|_i)2 T _ J g [y d)+<^od) ^_J_J (49c) 


where ^(1) and <p'o (l) can be found from equations 
(131) and (B4). Equations (44), evaluated at 
77 = 1 , become 

ro=i-rt +1 -a(i- v i+y« + ' (50a) 

e«=(D 0 +l)Ci(l-,I+ 1 )^ (50b) 

So=D 0 (D 0 + 1)^ 0 (1 — (50c) 

Solving for /+, C 0 , and rj b yields 


•4(^-0 


, /T77 W,, n ,Y, ns, 

+ Vi( a +<? " -1 ) + -«r 

rt Qu ( 

u A+iVW 

M. 1 S.V 


The pressure distribution is found from 




yUVH) 
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Tlie constants P 0 , Q 0 , and S 0 have the following 
values when <r=0,l. For <r= 0: 


Q tj- 


7—1 7 

_ 7+l 7-/3 

27-1 fi 


7+1 7-/3 


S 0= 


7 2 — 1 7 — /3 


[ (? 7— 5)- 


137-11 


(52a) 
(52b) 
sj (52c) 


For <r= 1 : 


Qu 


Po- 


7—1 7 


7+1 7— j8 


27—1 1 7(7—1) 

7+1 7— /?[+ (27— 1)(7+1)_ 


(53a) 

(53b) 


s '-2Fb)^[- (13T - nw ’ 

+2(j^^ g - (,-')(W10,+ |) ] (53c) 


Case 0=1. — For 0=1 it is known that 77^=0 so 
that an approximate solution is (from oqs. (50a) 
and (50b)) 

77^ = 0 — P 0 D 0 = (QJCo) 1 (54) 

This solution does not satisfy the boundary con- 
dition on <po(l). It gives more accurate results 
than do equations (49) to (51) evaluated for 
0 = 1 - 

Equations (54) correspond to the approximate 
analytical solutions of the zero-order blast-wave 
problem which are presented in references 3, 4, 
and 5. Equations (44) to (53) may be viewed 
as a generalization of the latter for 0<^1. An 
exact closed-form solution for the 0=1 case is 
presented in reference 10 (for cr=0, 1,2) and in 
reference 12 (for <r=2). 

Numerical results for .7=0,1 and various values 
of 0 and 7 are listed in tables I and II. 

APPROXIMATE ANALYTICAL SOLUTION FOR FIRST-ORDER 
PROBLEM 

The quantities y? 0 ), i/uM), an <I FJF Q will 

be considered as the dependent variables (similar 
to ref. 6). By using equations (40), the first-order 
equations (oqs. (16)) can then be written, re- 
spectively : 4 


Continuity: 

(— )' -(t)' +H--- +ii i)-- 

\y—<Po/ wo/ 00 \V — <Po Wo/ 


0 (55a) 


Momentum: 




+ 

I t\ 


Energy: 

- G [<-« jV* r ] " 

(55c) 


First, equations (55a) and (55c) will each be inte- 
grated as far as possible. 

Equations (55a) and (55c) can be written, re- 
spectively: 


(9-0 f/ 


\ 0 5 

V° +>; 

1 6l + & \ 

v 1 — <P0/ 


( 9-0 

7 1 

/ „ \/ 
0 7 ] 

V" PJ 

0(7 -/9+ 7 fi)/y \ 

fty 1 

v n—<Po/ 


Integration yields 

Integrating by parts and introducing the constants 
Ei and G x give the alternate expressions 



=-=+ 
V — <Po 

F t 

_ 7V>i _| 

Fo 

V Wo 


(l“h0)0o f™ ~ EM (56a) 

J “0 V <Po 


(56b) 


4 The operation [7(t-f0)~- 0] (eq. (55a)) — (1+0) (eq. (55c)) leads to the 
special integral 

0 ~-^ + (70+y-P) j-'-a+fl) - 5 =(coiist.)Pp 
q— <PQ ^0 *0 


where the constant can be evaluated In terms of the boundary conditions at 
3— I ■ 
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Substitution of equations (56) into equation ( 55 b) 
reduces the problem to the single dependent vari- 
able — <a>)- This is essentially what is done 
in appendix D to obtain the asymptotic form of 
the first-order solution near 0 o ^O. Since equa- 
tions (55) are equivalent to a linear third-order 
equation, three independent asymptotic solutions 
are obtained. When numerically integrating 
equations (16) (by the superposition procedure 
outlined in connection with eq. (19)), the asymp- 
totic solutions are required in order to proceed to 
7j b from a point near 77 * (since the equations are 
singular at 77*). 

An approximate solution of the first-order equa- 
tions (eqs. (16) to (18)) will be obtained herein 
by using the asymptotic solutions of appendix D 
and evaluating the arbitrary constants therein so 
as to give consistent values for the dependent vari- 
ables and their first derivatives at 17 = 1 . The 
asymptotic solutions which do not satisfy <^ 1 C 7 ;*) 
— 0 are neglected (in order to satisfy eq. (18)). The 
resulting solution will exactly satisfy the continuity 
and energy equations but does not satisfy the mo- 
mentum equation except at rj=l and 77 = 17 *. The 
cases 0 <^<^lj 0 = 0 , and 0=1 are treated sepa- 
rately. In each case the problem is reduced to 
a form requiring the simultaneous solution of four 
linear algebraic equations in four unknowns. 

Case O<0<1. — Let L u s= — pA u M it 1 = 

0(1— 0)<7 A u E lt2 = C u and Afi ti = hCi in appendix 
D. Equations (D4a) plus (D4b) can then be 
written 

=A l [- 0 +(Kl- 0 ) 9 e l >]+CMl +l, +( 2 -P)B l Ol 

y—<Po 

(57a) 


Define the coefficients of A x , B lt and C\ by 


&u = 0[(l— P)g~ 1] £>12 

&2i = l+2^gr £i 22 

hi = (y— P)+(2y— 0)0g b 32 

*41 =>3(1 — 0)9 &42 


I — A { (1 -\-2f3y6 0 ) + Ci0 s [l + (2 — j- >3) A^ 0 ] -)-3Bi0§ (57b) 

YO 

■^c = A 1 1 (7 — 0 ) + 0 (27 — 0 ) <j0 o ] 

+ CiA(27+70-0)*J + '+ (37-0)7?!^ (57c) 

where 

-1 M a>l -(7-0-1)^-^ 

9 1-0 Li.x 2(27-0)F o (i7*) 

i t 2 Pvl * 

n ~ E Xi2 — 2(1 +0) (27+07— 0)2^ (77 *) 

The term involving B v dl was added to equation 
(57a) so as to permit satisfying an additional 
boundary condition at 77 = 1 . The corresponding 
terms in equations (57b) and (57c) were found 
from equations (56) with 7ti = 6 ri = 0 . Four un- 
determined constants, A h B u C u and a u remain. 
These are evaluated so as to give consistent values 
for the dependent variables and their first deriva- 
tives a t 77 = 1 . 

At 77 = 1 , 0 O (1) = 1 and equations (57) give 
(recalling eqs. (17) and (B5)) 

^15 J)\±d \ = [— 0+0(1 — 0)</]/li + (2— 0)Bi-\-hC\ 

(58a) 

f>25 — &24®1= (1 + 20 g r )Al 1 + 3/?i+ [1 + (2 + 0 ) /*•](+ 

(58b) 

ho — h&i = [ (7- — 0 ) + (27 — 0 ) &g]Ai 

+ ( 37 — 0 ) 2 ?!+ ( 27+7 0— 0)ACi (58c) 

^45 — 644^1 — 0(1 — 0)<+h + 2(2 — 0)2?j + (1 -\- ff ) hC } 

(58d) 


2-0 

^13 = h 

3 

^23 — 1 + (2 + 0 )A 

3t-/3 

633 ^ (27+70— 0 )A 

2(2-13) 

^43 — (1 +0)^ 
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Equations (58a) to (58d) can then be written, 
respectively: 


611^11 + 612^1 + ^ 13^1 + 614^1 = ^15 
621^1 + 622^1 + 623+1 + 624^1 = 625 

631^1 + 632^1 + 633+1 + 634^1 = 635 


(59) 


641^1 + 642-^1 + 643+ 1 + 644^1 — 645 J 


which are four simultaneous equations for the four 
unknowns A u B u C u and a x . A knowledge of these 
four quantities, together with equations (57), com- 
pletely defines the first-order flow field. Numeri- 
cal results for < 7 = 0,1 and various values of p and 
7 are listed in tables ITT and IV. 


Case P = 0, — Let E 1A 


= C i} G lt2 = A l} and g = 


Mm 

Ei A 


= — ^ 1,s in appendix D. 
61 .2 

(D5b) become 


Equations (D5a) and 


V <Po 


(l+l A tl M+~ B 1 el 

J +| 7 Cigti+At (l -| 7 <l6l)+j BX 


Y (60) 


where 


K=—g 

7 t 3t 

* 31 = 1 -jV 

641 =— 2 <J 



&13 — 9 

, , , 3 

623 — 1+2 9 

, _3t 
O 33 — ~2 9 

b i3 =2g 


(62) 


When these quantities are used in equations (59), 
the unknowns A u B X) C u and a x can be found, thus 
defining the first-order flow. 

For <j=0, it can be shown that +li = (3 y+1)/2y, 
5i = 0, Ci = — 2 /( 7 — 1 ), and ai=l. Also, F x (t? & ) = 
(37+ 1 ) /[y (7 + 1 ) ] - These results are in exact 
agreement with those from an expansion (in terms 
of 1 /(MS) 2 ) of the oblique-shock relations for flow 
over a wedge. 

Numerical results are tabulated in tables III 
and IV. 

Case p-^ 1 1 — Let L Xi \= A\, E x 2 = +1 j Q — 
—Mi ti /Li t i f bndk = Mi t 2 lEi, 2 ^ appendix D. Equa- 
tions (D6a) and (Dfib) then suggest 

- 5£l -=/1,[- 1 + CP- 1) gen + M < +1 + Br 2 

V~ <Po 

^=^[l + |l-<r|(P+l)^] 

+ CA (h — y> /*®d) +pij^ 4 Bi6q+ 2 y 


— <j(y — 1)2 I/T 

^~3T(cr+l)+/[(7 + l)+oW] rt+I,/y 


The terms involving B X &1 were added so as to 
permit satisfying an additional boundary condition 
at 17=1. At 77 = 1 , equations (60) give (using 
eqs. (17) and (B5)) 


615 — 6 14 a 1 — — </+j+ J9 i + <7+ l 

6 2 o — 6 24 a l =— ^ f7^i+;^ B\ + ^1 +^ <7^+i 

&«-&hBi=( lJQ g y 4i+ 4 + P,+f gCi 

645 — 644(74= — 20 J 4 1 + 3 i? 1 + 2 (/+ l 


( 61 ) 


Define the coefficients of A u B u and +1 by 


§=+> {( 7 - 1 ) + [y(P+ 1 ) ■ - 1 ] gd % } 


+ ( 7 _ p+27 Uh Cl e p Q + 1 


7 P+ 37-1 

P+i 


whore 


27+< t -1 
7(0-+ 1) 


P^ +2 J 


(63) 


(T-IKt-1) 

( 47 Z 137 (17+8) Kg TU+f) 

g - 27 2 (<r+l)P[7(P + l)-l]Po(0) 

(Y-D(g-l) 

{3y+cry-2)PK 0 1 ' (<T + 1) 
! -27 2 ( < r+l)(P+l)(7P+27-l)Po(0) 


503S91 — 59 3 
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The terms involving Bfi^ 2 were’ inserted to permit satisfying an additional boundary condition. 
The solution then proceeds as in the previous cases and results in four simultaneous equations for 
the four unknowns A h B u C h and a x . Equations (59) apply with the coefficients of A u B u and 
Ci being given by 

l)ff— 1 

6ai=l + Il-H(^+l)ff 
bn=y-l+[y(P+l)-l]g 
b«=(P~i)Pg 

&I3— 1 

6*=(P+3)/(P+l) 

633= (tP+3t— 1)/(P+1) 

b 43 =P +2 


bi2 — h 

b 22 =l + [(P+2)h/P] 
h 2 =[(yP+2y-l)h]/P 
bi 2 =(P +1)A 


r (64) 


Equations (59) can be solved for A X} B x , C X) and 
Ci. Numerical results are given in tables III 
and IV. 

NUMERICAL RESULTS AND DISCUSSION 

The zero -order and first-order problems were 
solved both by numerical integrations of the 
equations of motion and by the approximate 
method. The results are tabulated in tables I 
to IV for cr= 0, 1 and various values of 0 and 7. 
The results of the numerical integrations of refer- 
ence 10 are included in tables IT and IV. The 
quantities F 0 (n&)> FiiVb), ??&> and a x are used in 
equations (21), (22), (32), and (33) to find pressure 
distributions and shock shapes for the class of 
bodies considered in the present report. 

With regard to the zero-order problem, tables I 
and II show that the approximate method is in 
good agreement with the numerical integrations 
for 0 near zero. As 0 approaches one, the ap- 
proximate method becomes less accurate. At 
0— 1, however, the approximate solution is again 
accurate since the appropriate value >7*= 0 is auto- 
matically imposed and there are only two free 
constants, as opposed to the three free con- 
stants in the 0^ \ cases. In general, the approxi- 
mate method is accurate when the shock is rela- 
tively close to the body (i.c., is near one) 
so that, for a given 0 , the approximate solution 
is most accurate for values of 7 near one. The 


estimates for rj b tend to be more accurate than 
those for F 0 (r/ & ). The variation with 77 of the 
dependent variables is plotted in figure 2 for 
or=l, 7—1.4. Figure 2 is based on the approxi- 
mate zero-order results. Corresponding figures, 
from an exact integration of the zero-order equa- 
tions, are presented in reference 10. 

The accuracy of the approximate solution of 
the first-order equations can be deduced from 
tables III and IV. Again, the approximate solu- 
tion tends to be more accurate for 0 near zero and 
for 7 near one. The accuracy of the first-order 
approximate solution is less critical than that for 
the zero-order flow since the former is a perturba- 
tion quantity. Thus, if the first-order solution 
represents a 10-percent correction to the zero-order 
flow, and, if the approximate first-order solution 
is 10 percent in error, the latter would represent 
only a 1 -percent error in the over-all flow. 

The numerical integrations of the first-order 
problem, which are reported in reference 10, ap- 
pear to be in error, particularly with regard to a x 
(see table IV). Note that for a~ 1 , /3= 0, 7=1.4 
the present approximate method and numerical 
integrations both give ax^O.48, whereas reference 
10 gives ai~0.40. These can be compared with 
the value aj—0.47 from the cone results of refer- 
ence 2 (see table IV), indicating better agreement 
with the present results than with reference 10. 
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GENERAL CHARACTERISTICS OF FLOW FIELDS 
ASSOCIATED WITH POWER LAW SHOCKS 


It was previously shown that, in the limit 
1/(M$) 2 - >0, a shock shape of the form E 0 =x m 
reduces the hypersonic-slender-body equations to 
a set of ordinary differential equations with 77= r/7? 0 
as the independent variable. The alternate shock 

shape parameter 0=|j> 0+1) was also 

introduced. Emphasis was placed on the range 
O<0<1. In the present section, the general char- 
acteristics of the zero-order flow fields associated 
with different values of 0 (including 0<O, 0>1) 
will be discussed. 

With cr=0, equation (39) becomes 


i , (y+ir^yel+e / fl„\ 
aj_ _4 (7 = l) 7 W 7 V yB'J 
(e'o ) 2 7 . ( 7 +l) 7+l yOl+* Op 

2y(y — l) 7 (d' o y V d' 


(65) 


*o(l) = l flo(l) = (7+l)/(Y-l) J 


Equation (G5) is singular when 0 o =O. The cor- 
responding value of 7 7 is denoted by It can then 

be shown that for 0 O — >0, 77— >17$, equation (65) has 
the following asymptotic forms: 


0<7* 0o~ K 0 (v— VbV 0 | 

0=7: 0 O 

y 

7<^<27: 0 O ^ Kq/(— 7i) 0 ~ y 

where A" 0 and L Q are positive constants. Asymp- 
totic solutions were not found for 0>2y. Also, 
recall that 77*= 0 for 0=1. 

I.'' 




V 

Figure 3. — Variation of d 0 with tj for <r=Q } 7 = 1 . 4 , and 
various values of j3. 


°f 

8 - 

6 - 


Curve 

0 

b 

c 

d 

e 

f 


— r- 1-' 


0<O 

0-0 

O<(0 = 3/4) < I 
0 = * 

l<{0 = 1 . 2 ) < 1 .4 
1.4 < 0< 2.8 

Numerical integrations 
of zero-order equations 
Schematic 


-.6 



-2 


11*=— 1 


Vb=— 00 J 


r ( 66 ) 



(a) 


r 



<c ) 




(a) 0 <O. 

(b) /3— 0. 
(c) 0</3< 1 . 

W) 0— I . 
(e) 1<0<7- 
(0 7</3<2 7 . 


Figure 4. — Schematic representation of flow fields for 
o '— 0 and various values of 0 . 


The formal solutions of equation (65), for con- 
stant 7, can now be represented in the 0 O ,7? plane 
for various values of /I This is done in figure 3. 
The corresponding physical flows are given in 
figure 4. The value of 0 defines the shock shape. 
The corresponding value of r} b defines the body 
shape (since r b —’q b B {) ), For rj b positive, the body 
is in the first quadrant of the x,r plane (figs. 4(a) 
to (c)), while, for 7 ] b a finite negative number, the 
hotly is in the fourth quadrant (fig. 4(e)). The 
body and shock arc similar in shape except for 
7/^=0, — 00. The drag associated with these 
flows can be found from equation (29). For 
0>1 it is seen that the drag is infinite at x — 0 
and thus cannot correspond to a physical situation. 

Except for the 0=0 case, each flow field in 
figure 4 is not everywhere consistent with the 
original use of the hypersonic-slender-body equa- 
tions. In particular, conditions at #=0 violate 
the hypersonic-slender-body approximations for 
all except the 0=0 case. However, in each case, 
there is a region in the vicinity of the shock 
(excluding x >0) for which the flow field is con- 
sistent with the approximations. If a streamline 
in this region is taken to define a body shape (see 
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sketch (c)), then t lie flow external to this stream- 
line is accurately described, at least initially, by 
the solutions obtained in the body of the report. 
Thus, even the cases yield flows, portions of 
which are physically realistic and consistent with 
the hypersonic-slender-hody approximations. 
The streamlines downstream of the shock are 
defined, parametrically (with 77 as parameter), by 


x 


(£>" 



(67) 


where ( x s ,r s ) are the coordinates of a streamline at 
its intersection with the shock. The ratios x/x s 
and r/r s are thus functions only of 77 (for specified 
c r, 7, and /3). Note that 0 O and i] are not constant 
along streamlines which do not pass through the 
origin. The downstream extent to which such 
solutions are valid has not been resolved. 



A similar discussion can he made for the <r~ 1 
case. Considering equation (39), the asymptotic 
form for as 0 o -->O, r}-->rj bj is, for /3<1, 

7 

8 0 ~K 0 (i] 2 —niy~ e ( 68 ) 

where 77^=0 for 0=1. The asymptotic form 
when has not been determined. Since the 
flow is axisyminetrio, 77 cannot be negative and the 
integral curve is confined to the first quadrant of 
the plane. There is a singularity at 77=0. 
If a streamline, other than 0 o =O, is taken to define 
a body, then the external flow corresponds to the 


flow about an open-nosed body of revolution. 
The leading-edge angle has a finite nonzero value 
for such flows. 

CONCLUDING REMARKS 

An approximate analyt ical method has been pre- 
sented for obtaining the zero-order and first-order 
solutions for hypersonic flow over slender blunt- 
nosed bodies following a power law variation. 
Wedge, cone, and constant-energy flows are in- 
cluded as special cases. The solutions arc found 
within the framework of hypersonic-slender-body 
theorjL 

The approximate solutions are compared with 
numerical integrations of the zero- and first-order 
problems. The agreement is generally good, par- 
ticularly for 0 near zero and 7 near one. The 
shock is relatively close to the body for the latter 
cases. Sufficient numerical results have been 
tabulated to permit estimates of the accuracy of 
the approximate method for various combinations 
of and 7. 

The general characteristics of flow fields asso- 
ciated with power law shocks have also been dis- 
cussed. It is pointed out that values of the shock 
shape parameter j3 outside the range O<0<1 also 
give rise to flow fields, portions of which are 
physically realistic and consistent with the 
hypersonic-slender-body approximations. 

The accuracy of the approximate solutions 
could be improved by following more complex pro- 
cedures. For example, the flow near the shock 
can be expanded in a Taylor’s series about 77=1, 
and the flow in the vicinity of the body can be 
expressed in terms of the asymptotic solutions, 
each multiplied by an arbitrary constant. These 
constants can then be evaluated by matching the 
Taylor’s series and the asymptotic expressions at 
some point between the shock and body. The 
Taylor’s series then represents the flow between 
the match point and the shock, while the asymp- 
totic expressions represent the flow between the 
body and the match point. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, November 17, 1958 
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APPENDIX A 


SYMBOLS 


A u B u O 

constants (eqs. (57), el seq.) 

y 


constant (eqs. (11)) 

5 

K 

constants (eqs. (17), (B5), (58), et 



seq.) 

8 b 

r D 

drag coefficient, (eq. (25)) 


c. 

pressure coefficient (eq. (20)) 

€ 

C<); D 0 

constants (eqs. (44)) 

V 


specific heat at constant volume 


D(I) 

forebody drag 

00 0?) 

F v W, Fi(n) 

pressure similarity variable (eqs. 



(11)) 

1 

T 

integral (eq. (30)) 

p 

K 0 

constant (eq. (42)) 

a 

L 

M 

streamwise length of body 
free-stroam Mach number 

<PiM, <pi(v) 

m 

shock shape power law exponent 

Pi O j Qi), Sq 

(eq. (6)) 

constants (eqs. (49)) 

tAoWi $ i(v) 

V 

pressure 


<1 

dynamic pressure, pJiP ^/ 2 

Subscripts: 

lip) 

lateral displacement of shock 

b 

ffoW 

lateral displacement of shock in 

s 


limit as 1/ (A/5) 2 — >0 

0 

r b (7) 

body shape 

1 

T 

temperature 

oo 

u,v 

velocities in f x,r ) direction, re- 

Superscripts 

n 


spectively 


streamwise and lateral coordinates, 


respectively 

( ) 

P 

alternate shock shape parameter, 



Bfi-OlA'+i) 

( y 


ratio of specific licats 
characteristic slope, equals 1 ? 0 (L)/L 
for flows with power law shocks 
characteristic bod}" slope, T b (L)jL 
perturbation parameter, 
lateral coordinate similarity varia- 
ble, rj 7? u (eqs. (8)) 
zero-order stream function simi- 
larity' variable (eq. (3G)) 
x (eqs. (8)) 
density 

0,1,2, for planar, cylindrical, and 
spherical flow, respectively 
lateral velocity similarity variable 
(eqs. (11)) 

density similarity variable (eqs. 

(ID) 

quantity evaluated at body surface 
quantity evaluated at shock 
zero-order solution 
first-order solution 
undisturbed frec-stream value 

barred quantities are dimensional 
unbarred quantities are nondimen- 
siorial (eqs. (1)) 

primes indicate differentiation with 
respect to rj 


The approximate analytic solutions require a 
knowledge of derivatives of the dependent vari- 
ables at tj= 1. These are summarized as follows. 

From equations (13) and (14) it can he shown 
that 

^e)=(^qhy2 [ ; s(t+i)(*+i)/3-4<t 7] (Bi) 


APPENDIX B 

DERIVATIVES OF DEPENDENT VARIABLES AT v 

2 


K0) = 


"(7-1) (7+ H 2 


[(27-l)(7+l)(»+D/S 




(7—1) 


— 2<ry(7 — 1)] (B2) 

,[3(7+I)(T+l)/9-2ff(7-l)l (B3) 
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+ 




P+ 


ry(y—l),y— 


7+1 2 

Similarly, from equations (1C) ami (17), 


-]n (i) -[(2-Y)( ? ~-)^ 


<t(y — l) 2 . 7 + 1 


7 + 1 


] ,o0 ) 


iTnT (- —)' = ( +i u +i )» {\ y iT K0 } (i) +7<r ' fT+ 1} (ff+ 1 ^ 1 (7+ 1)2 
0<j(l) \V — iPoA-i (<r+l)(7+*l V.L ^ 


7—1 


, (7+1)11 ,.(!)_ 

' r 2(7-l)J 1-^0) 


+7^(1)+^ 

Lt— i y— 


2a ~|/ ?T 1 ( 1 ) 7 — lpy + l 

'0)“ 2 Lt -1 


**(!)- 


(7 + l)(<r+l)^ n ^(l) \ _ 7 y 
(t— i) 2 JMVJ 45 44 1 


(B4) 


(Bo) 


APPENDIX C 


BOUNDARY CONDITIONS AT 17= 1 FOR FIRST-ORDER PROBLEM 


In the first-order problem, tin* shock is located 
at R = £™(l + tfi€) or, equivalently, at r)=l + aie 
(sec sketch (cl)). It is desired to specify boundary 
conditions at 77 — 1 . 



Let Q be any quantity whose value is known at 
the downstream side of the shock (designated by 
subscript s ) . Then, expanding in a Taylor's series 
about the shock at constant x yields 

(«),-!= (®.+(^) f (ffo-ff)+ • • • 

- ( ®-'(f)+ ,+ '" (C,) 


Consider Q to be the vertical velocity. Equation 
(Cl) becomes 


i\=i =v. 


8 (>»)."“ 


(C2) 


But, from equation (4e), v,=^j +/ or 

.^+ir-[. + .( 2 -+»,-i)] (C3) 

since tl 7?/d j c = |]+e( ^ ? ”) «,]• Equation 

(C2) then becomes, by letting isfflf' 1 (v , u+«^j) 


«>(!) + «w(l) = 


7+1 


-±Y1 

|_ \ m in vj 


so that 

*<»=4i {“i + “( 2 


'V?J(l)cfi€ 

(C4) 


1 — m 7 r 


m 


- ,0(1)] \ 


(C5) 


Expressions for 7*1(1) and ^i(l) are found similarly. 
The results are given in equations (17). 
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APPENDIX D 

ASYMPTOTIC SOLUTIONS OF FIRST-OROER EQUATIONS 

For A T =0, P= l } E ltl = G itl ^ Q: 


Asymptotic solutions of the first-order equa- 
tions, near 0 o asO, are found herein. 

From equations (56), the continuity arid energy 
equations show 


*P\ 


V V^o J V <Po 


F_ <p i 

F 0 ~ y v~ 


<Po 


•( 7 +7 v-m f 


On (1+ ^ ) 


<PI 


y—<po 


Consider solutions of the form 


<p i 


V ~<Po 


-9q d - • ■ •) 


(Din) 

■ d0 o -j-G+0o 

(Dtb) 


(D2a) 


For the corresponding values of 

\pi and F x are 


*1 »wV-l 






p,+ 


V+P+l 


F, 

F 




(3 ¥+P-3 

~3 


!/,«?+... )+/5i*8 
(D2b) 


7(A r +P)+7— /3 , , p 
+ N+P-p' lUl9 °' 


'■] 


f<W8 (D2 c) 


The appropriate values of A 7 and P can be found 
by substituting equations (D2) into the momen- 
tum equation (cq. (55b)) and considering 0 O « 0. 
Tiiis will be done for £=(), and fi=l, 

respectively. 

CASE O<0<J 

In the vicinity of #y«0, the momentum equation 
becomes 

(D3) 

Substitution of equations (D2) into (D3) yields 
the following three independent solutions: 


<p i, 


V~<Pq 


- Pi.i+M i( i0 o + • ■ • 


^1,1 PiA t 2il/j ] 

d + 1-/° 


^0 

Fi,!_ 7 0/ 27-3 

T" 1 ' 1+ ‘ 1 -p /l>10+ 

Cy— g— 1) (1— 

Pi. i 2(27— P)F 0 (vt l ) 

For .¥=3, P=l, ^1.2—0: 


r (F4 ii) 


.2 

»7 — 


^o(b+^^l,2^0+ * • 0 


— 0o[P^2+ (2+/?)d/ 1 o0 o + . . ,] 

ro 

“tP— 0o[Op $)d/ x oftid* . . .] 

P lt3 “2(l+^)(2 7 +^ 7 -/3)Po(^) 

For tf=- (T-fl/ 7 , P=l, ^ 1 , 3 -^, 3 = 0 : 

’j- W) 


(D4b) 


^ 1, 3 
1^0 

F.3 

Fo 

•U,3 


a ~ Z -> g f ZZ^Lhl 1 !/■ a I 

00 Lw3-/3 + 0(i-7r /,,,fl#+ ' ' + 


=00 


t 


0+ jfc-7) * 0+ • ‘ ’] 


(y-\)p 2 v\^ 


L u 3 27(7+73-3) F(W 


(D4c) 

For solutions satisfying v 5 ! ( 77 O = 0 it is necessary 
to omit equations (D4c). 
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CASE )S = 0 

The momentum equation can be written 


When the boundary condition — 0 is to be 

satisfied, equations (D5e) are omitted. 


— ((7+1) 


00 d[<Pi/(v— <Po)] , v „ d (F } /F 0 ) 
do ddo TT 


-jp (l+ 2 ff) 



where 

9'**{v + \)rfK 0 ~(c+ 1),' (£t[) [+' n(^)]' /r 
The following three independent solutions result: 


For A T = 0, A=2, (?i il =0: 
-^=0+M 1>1 flg+ . . . 

V — <Pi} 

+=A,.+ A/ m 81+ . . 

ro " 

^=0 + y A/,.,flg+ . • • 


r CD 5a) 


M,.,_ <r(Y— 1)2 1/7 

A,., 37(<7+l)^[(7-t l)A 0 („ 6 )]"+n/r J 


case 0=1 


For this ease, jj ~ ^*=0, - , u 


7-1 Oo 7<r+I 


0'y 7(<r+l) (O’o) 2 7(^+1) 

2y + (j~ 1 (7-0 (g~0 2^ -}-g — 1 


and 


■ A’ui? 


T-l 


A 0 


70+1) 


^0 


70+1) 


Recall that 



I he momentum equation can he written 

7 1 „ Wi)l | jr d(A,/A„) 

M d^ +F ° 


( 7 — l)(ff7— 7+4 ) 
27 2 (c+1) ’ 


3t+T(7— 2 

2t 2 ((7+ 1 ) 




The following three independent solutions are 
obtained : 


For A 7 =0, P=2, /? 1|2 = 0: 

J?L? -=o+ii/ lll flg+ . . - 

V~<Po 

+=o+|m,, 2 ^+ . . . 

ro 2 , 

A< . 2 ri I ^ 1 / n 2 I 

”^T" — 'ri,2 + "2* dli |2 0 O + 




M,, 2 = o~(7 — 1)2 1/Y Mm 

C Ii2 37(<r+l)^[(7+l)F 0 (W] (v+,,/t_ Am „ 

(Dob) 

For A 7 = — 1, A= 2, C 1 , 3 =A j , 3=0: 


J ^=0 U - 1 (A,3 + M 1 ,30§+ • • •) 

7 7~ «A) 

^^==0(7 1 (0+2il/ i3 0 2 + . . .) 

ro 

X 

(O+27A/1 3 0 2 + . . .) 

1 0 

M,.,_ < 7 ( 7 — 1)2 ,/T _ 3 A/ m 

/ 1I .," _ 27(a+l)^'[(7+l)Ao(,*)]"+‘>/’- 2 Am 7 

(D5e) 


For A T =(), A=(2 7 +<7 — l)/[ 7 (<r+l)]: 

( 1 ) < 7 = 0,2 (i.P., 7 Vl): 


Vl.I 

V—<Pu 


=n. I +(p-i)A/,x+--- 


Pu 

Pa 

Am 

A n 


+=-A I . 1 +(A+l)A/ J ,,C+... 

Pa 


= -(7-1)L,. i +[ 7(P+1)-1]M i . 


P 

1^0 " 


(4 7 2 — 1 37— <77+ 8) A 0 


(7-1) o-l) 
70+1) 


Ai,i 2 7 2 (<7+1)A[7(A+1)-1 ]A u ( 0) 
(2) <7=1 (i.e, A=l): 


77—^0 


o. li ,+o+ . . . 


+— r,,+o+... 

ro 

~y = — (7 1 ) Am+(?i_ < 00 + ■ • • 
<0 

Am V-7t+4 
Am 2 7 2 F o (0) 



24 


TECHNICAL REPORT R— 1 5~ - NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


Those can bo pul in the unified form, valid for 

(7 = 0, 1,2: 

• •• 

17 — V»o 


!^I = ._7: l . 1 +|l-(r[(P+l)M l X+ 

ro 

"o 


y (ld6a) 


-1PZX+--- 


(7-1) (°-i) 

d/, .i U7 1 2 — 1 dy — <77 -f 8 ) K q l (a+1) 

Z,., _ 2 T 2 (ff+l)P[7iP+ l)"-1]Fo(0) 

✓ 


For A T =1, P=( 27 +<r-l)/[ T (< 7 +l)], £,.*=0: 


^_=fl 0 (0+M,.A P + • • •) 


v u V 

y—<Po 

- i ^ 1 , U 1 ’ ' / 


^, 2 + p + 2 * 5 6 a/,a p + • • •) 

■t 0 

^ ()+ 7 /, + 2y ZL 1 MiJ r + . . 

M, , 

(3 T + (77 — 2 ) P (Zlo) 7(o+1) 

to 

ii 

to I 

; (< 7 +l) (P+ 1 ) ( 7 P + 27 — l)Po( 0 ) 


(Deb) 


For A T = — ( 7 — 1)/t, P=( 27 +o— l)/[ 7 (( 7 +])], 
Ei,i~ fd, a" 0 : 


7-1 


-^=0 U 7 (Z M +A/ lt A p + • • •) 

y—vo 


Po 


! =0U 


/ Z, 3 

Vl-2 


7P-) 

7P+I- 


-27 


A/ M 



7p =<?u 7 ( 0 +^iZ^ m *> 3 ^+ • • • 

C7 — D(<r — 1) 

Mi 3 r(37+^-2) + (2 7 -l)(y-l) (<ry-y+4) ( 7 f > ~M~ 27 )/vo + j 

L^r L ' 2 7 2 (<z -p 1 j (2 7 — 1 ) ‘ V t / J yP(yP+ i-y)PM 


r ( D6c ) 


When the boundary condition <pi0? 6 )=0 is to be satisfied, equations (D 6 c) arc omitted. 
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